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Design  of  a  Multivariable  Helicopter  Flight  Control  System  for 
Handling  Qualities  Enhancement 
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New  handling  qualities  specifications  are  currently  being  developed  tor  attack  helicopters.  Most  unaugmcnled 
helicopters  will  not  meet  these  specifications  and  feedback  control  is  necessary  to  improve  handling  qualities  so 
that  safe  operation  close  to  the  earth  in  poor  weather  conditions  and  or  at  night  is  possible.  In  this  paper  a 
methodology  for  the  direct  design  of  helicopter  flight  control  systems  which  meet  handling  qualities  specifications 
is  presented.  This  methodology  uses  full  state  feedback  to  place  closed  loop  eigenvalues  to  achieve  bandwidth 
specifications  and  to  shape  closed  loop  eigenvectors  to  decouple  lateral  and  longitudinal  responses  to  control 
inputs,  full  state  feedback  requires  that  all  state  variables  he  known:  however,  only  angular  rates  and  normal 
acceleration  are  measured  In  sensors.  Thus,  a  state  estimator  is  required  in  the  feedback  loop  in  order  to  convert 
sensor  outputs  to  control  inputs.  This  estimator  is  designed  using  eigenstructure  assignment  so  as  to  achieve 
loop  transfer  recovers.  Design  of  a  feedback  system  for  use  in  precise  hovering  control  for  a  modern  attack 
helicopter  is  used  to  illustrate  the  method.  Control  law  synthesis  is  accomplished  using  an  eighth  order  model 
which  includes  only  rigid  body  modes.  Control  law  performance  is  evaluated  using  a  V7th  order  model  which 
includes  rigid  bodv .  actuator,  rotor,  sensor,  and  flexure  dynamics.  It  is  found  that  a  notch  filter  must  he  added 
to  the  design  in  order  to  eliminate  a  high  frequency  instability.  Once  this  is  accomplished,  both  the  time  and 
frequency  response  characteristics  of  the  augmented  helicopter  are  much  improved  compared  with  the  unaug¬ 
mented  helicopter. 


Notation 

Scalars 

/  =  \  -  1 

/  =  arbitrary  number  of  eigenvalues 

in  =  number  id  controls 

»  =  number  ot  states 

ii  =  normal  acceleration  of  e  g. .  g’s 
p  -  roll  rate,  rad  sec 

q  =  pitch  rate,  rad  sec 

i  =  yaw  rate,  rad  scc 

v  -  Laplace  operator 

n  =  forward  velocity,  ft  see 

=  lateral  velocity,  ft  sec 
a  ~  downward  velocity .  It  sec 

:  =  i-th  transmission  zero 

u ,  -  collective  pitch,  deg 

to  =  longitudinal  cyclic  pitch,  deg 


This  is  a  revised  version  of  a  paper  presented  at  the  4ird  Annual  Forum 
of  the  American  Helicopter  Society .  Si  Louis.  Mo  .  May.  14X7 


tix  ~  lateral  cyclic  pitch,  deg 
bj  =  tail  rotor  collective  pitch,  deg 
«r(/4>  =  minimum  singular  value  of  matrix  A 

<f(A)  =  maximum  singular  value  ol  matrix  A 

f)  =  pitch  angle,  rad 

q  =  roll  angle,  rad 

Vectors 

u  -  control  vector.  |u,.  to.  to.  i<4|' 

it,  =  command  control  vector  trom  control  nuxer 

t  =  state  vector.  |u.  r.  a  .  p.  q.  r.  if.  f)|' 

x  =  estimate  of  state  vector 

v  =  measurement  vector.  | p.  q.  r.  n_|7 

=  pilot  input  command  vector,  (collective,  longitudinal, 
lateral,  directional]' 

Matrices 

A  =  open  loop  dynamics  matrix 

B  =  control  distribution  matrix 

C  =  measurement  distribution  matrix,  state  vector 

D  -  measurement  distribution  matrix,  control  vector 
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£(v)  =  multiplicative  error  matrix.  GY  *(.v>  G.-m  -  / 

<7(0  =  open  loop  transfer  matrix.  C't/.v  -  A)  'ft 
fti.o  =  full  state  loop  transfer  matrix.  A'i/s  -  A )  'ft 
H  =  eontrol  mixer  for  pilot  commands 
A'(.i)  =  compensator  transfer  matrix.  A‘(/.v  -  A  *  ftA  -t 
AC  i  'L 

A'  =  feedback  gain  matrix 

L  =  estimator  gam  matrix 

Superscripts 
T  =  transposed 

-  I  =  inverse 

Subscripts 

K  -  8th  order  model 

37  =  37th  order  model 

Introduction 

ecently  a  number  ol  papers  have  appeared  which  discuss 
the  application  ol  various  modern  feedback  control  design 
techniques  to  helicopter  flight  control  synthesis.  In  the  past, 
classical  single-input-single-outnut  (SlSOt  frequence  response 
techniques  have  been  used  to  design  control  laws  for  helicopter 
llight  eontrol  systems.  However,  since  helicopter  responses  to 
control  inputs  arc  highly  coupled,  helicopter  dynamics  are  char¬ 
acterized  by  multi-input-multi-output  (MIMOl  mathematical 
models  and  use  ol  classical  SISO  techniques  max  require  a 
great  deal  of  time,  consuming  trial  and  error  effort  Modem 
Ml  MO  techniques  are  well  suited  to  the  design  ol  control  laws 
tor  helicopters,  and  numerous  papers  have  appeared  which 
describe  such  applications  These  include  linear  quadratic  reg¬ 
ulator  theory  (Rets.  I  -  2 ».  multivariable  loop  shaping  i Ret.  3). 
model  following  i  Ret  4i.  optimal  output  feedback  design  (Ret 
5t.  and  H‘  techniques  t Ref.  b).  The  advantages  and  disadxan 
tagcs  ot  some  of  these  techniques  are  discussed  in  Ret  5  The 
objective  of  this  paper  is  to  present  the  application  of  another 
technique,  eigenstructurc  assignment,  to  the  design  of  a  heli¬ 
copter  flight  control  system 

Eigenstructurc  assignment  is  a  technique  for  synthesis  of 
feedback  control  laws  which  allows  the  designer  to  directly 
place  closed  loop  eigenvalues  and  eigenvectors  in  specified 
configurations  These  configurations  are  selected  so  that  the 
closed  loop  response  characteristics  of  the  controlled  helicopter 
satisfy  handling  qualities  specifications.  In  this  paper  eigcn- 
structure  assignment  is  used  to  synthesize  control  law  s  for  pre¬ 
cise  control  of  forward,  side  and  vertical  velocity  and  yaw  rate 
lor  an  attack  helicopter  in  hover  It  is  desired  to  have  pilot 
longitudinal  stick  commands  correspond  to  forward  velocity . 
lateral  stick  to  side  velocity .  collective  to  vertical  velocity,  and 
pedal  position  to  yaw  rate.  Coupling  between  longitudinal, 
lateral,  vertical  velocity  and  yaw  modes  is  to  be  minimized. 
The  closed  loop  bandw  idth  must  be  large  enough  to  insure  crisp 
response  to  pilot  inputs  and  the  closed  loop  dynamic  response 
should  be  stable  in  the  presence  of  errors  in  the  design  model, 
due  to  effects  such  as  unmodeled  dynamics,  nonlmcaiities.  and 
variations  in  parameters 

The  rest  of  the  paper  is  divided  into  five  parts  First,  the 
mathematical  models  used  for  both  control  system  design  and 
evaluation  are  discussed  Next,  the  performance  requirements 
which  the  closed  loop  helicopter  must  meet  arc  described  The 
design  of  the  control  system  is  outlined  in  the  third  major 
section  and  the  evaluation  and  modification  of  the  design  is 
given  in  the  fourth  section.  The  last  section  consists  ol  con¬ 
clusions  and  suggestions  for  additional  work. 

Mathematical  Models 

The  helicopter  modeled  in  this  study  is  a  modem  attack 
helicopter  similar  to  the  YAH-64  (Ref  3).  The  control  laws 
are  designed  using  an  eighth-order  rigid  body  model  and  eval¬ 


uated  using  a  thirty -seventh-order  model  w  hich  includes  ac¬ 
tuator.  rotor,  sensor,  and  flexure  dynamics  Mam  rotor  collective 
pitch,  longitudinal  cyclic  pitch,  lateral  cyclic  pitch,  and  tail 
rotor  collective  pitch  are  the  control  inputs  Three  body  rate 
gyroscopes  and  an  accelerometer  which  measures  normal  ac¬ 
celeration  are  used  as  sensors.  The  mathematical  model  used 
is  semi-empirical  and  was  developed  for  the  hover  llight  con¬ 
dition.  The  linearized  rigid  body  equations  ot  motion  are  ex¬ 
pressed  in  standard  state  variable  form  as 

,i  =  ,4.i  -  Bit  III 

y  =  r.i  a-  Du  ( 2 1 

The  A.  ft.  C.  and  D  matrices  were  obtained  by  numerical 
linearization  of  a  nonlinear  analytical  model  ol  the  helicopter 
and  are  given  in  Table  I  The  actuator,  rotor  and  sensor  flexural 
dynamics  models  used  are  given  in  Fig  I  and  were  obtained 
from  flight  test  data.  Sensor  outputs  were  measured  and  re¬ 
corded  for  pilot  inputs  at  various  frequencies.  The  transfer 
functions  shoxxn  in  Fig.  I  were  obtained  by  fitting  assumed 
forms  of  the  input  and  output  transfer  functions  coupled  with 
the  model  given  in  Eqs.  ( 1 1  and  1 2 1  to  the  measured  data  and 
numerically  adiusting  time  constants,  damping  factors,  and  nat¬ 
ural  frequencies  until  the  frequency  response  ol  the  model 
matched  that  of  the  helicopter.  The  \ alidity  of  this  type  ol 
model  for  flight  control  design  was  verified  in  Ret  3  in  which 
a  flight  control  system  was  designed  using  a  model  similar  to 
that  given  in  Table  1  and  Fig  I  This  control  system  was 
successlully  flight  tested. 

The  open  loop  eigenvalues  and  non-dimensional  eigenvectors 
of  the  design  model  are  given  in  Table  2.  Non-dimcnsionali 
zation  of  the  state  vector  was  achieved  by  dividing  the  linear 
velocities  by  25  It  sec.  the  angular  rates  by  20  deg  sec.  and 
the  angular  displacements  by  20  deg  These  are  the  maximum 
values  of  the  state  variables  expected  during  hover  maneuvers 
The  control  variables  were  non-dimensionalized  by  dividing 
by  their  maximum  values  ot  4  deg  collective.  15  deg  longi¬ 
tudinal  cyclic.  8.75  deg  lateral  cyclic,  and  18.5  deg  tail  rotor 
collective  pitch. 

Examination  of  the  eigenvectors  and  eigenvalues  in  I  able  2 
indicates  significant  coupling  between  lateral  and  longitudinal 
modes  Both  forward  and  side  velocities  exhibit  low  frequency 
instabilities  and  as  is  discussed  m  the  next  section,  the  band¬ 
width  in  the  pitch,  vertical  velocity .  and  yaw  directions  is  not 
large  enough  to  guarantee  level  I  handling  qualities  In  addition 
to  the  modal  coupling,  examination  ot  the  control  distribution 
matrix  ft  in  Tabic  I  reveals  strong  control  coupling  between 
lateral  cyclic,  longitudinal  cyclic,  and  tail  rotor  collective  Only 
main  rotor  collective  is  relatively  uncoupled  Time  histories  of 
open  loop  responses  to  a  lateral  cyclic  step  input  are  shown  in 
Fig.  2  The  response  of  the  unaugmented  helicopter  is  indicated 
by  solid  lines.  It  can  be  seen  that  the  responses  about  all  axes 
become  very  large  due  to  the  unstable  eigenvalues  and  signif¬ 
icant  coupling  which  exists  between  modes. 

The  dynamic  responses  for  the  helicopter  modeled  in  this 
study  are  typical  of  most  high  performance  helicopters  Sim¬ 
ulations  and  flight  tests  have  shown  that  even  experienced  hel¬ 
icopter  pilots  arc  unable  to  accomplish  relatively  simple  hover 
tasks  in  conditions  of  degraded  visual  cueing  and  or  divided 
attentions  tasks  with  such  typical  helicopter  dynamics  (Refs. 
7.  8).  However  flight  tests  with  variable  stability  helicopters 
have  shown  that  stability  augmentation  is  an  effective  method 
for  compensating  for  missing  v  isual  cues  and  for  use  in  situ¬ 
ations  in  which  the  pilot  must  devote  a  significant  amount  of 
time  to  tasks  other  than  piloting 

Performance  Requirements 

Three  response  types  have  been  proposed  to  quantify  mission 
oriented  rotorcralt  handling  quality  requirements  (Refs.  7-4). 
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Table  I  Stale  and  (lain  Matrices 


-  .0286 

- .0637 

.0205 

.2240 

7.9700 

-  .2570 

.0000 

—  32.0000 

.0774 

-.2310 

.18)54 

-8.2400 

-  1.0300 

-  1.6400 

32.0000 

.1640 

.0046 

-  .0257 

-.2610 

-.3740 

2.2500 

2.1900 

1.6000 

-  3. 2818 1 

.0074 

— .0500 

.0095 

-  2.70000 

-.1340 

-  .6620 

.0000 

00(H) 

.0047 

.0118 

.0002 

- .0042 

- .7500 

.0244 

.0000 

.000(1 

.0034 

-.0044 

.0008 

-  1.0500 

.4130 

-  .4000 

.0000 

.0000 

.0000 

.0000 

.0000 

1.0000 

— .0051 

.1030 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.9990 

.0494 

.0000 

.0000 

.4350 

.5760 

-.1140 

-  .0004 

-.1580 

.1360 

.4410 

.2820 

-4.27(H) 

.0575 

- .0250 

.0012 

- .0438 

-  .0600 

.6470 

.0800 

.0072 

-.1010 

-  .0900 

-  .1)019 

.0800 

.0047 

.2000 

-  .0455 

.0000 

.0000 

.0000 

.0000 

.0000 

.IKItHI 

.0000 

.(loot) 

l  = 


.0000 

.1888) 

.1888) 

1.(888) 

.0000 

.0000 

.00(H) 

.0000 

.0000 

.18)00 

.0000 

.0000 

1 .0000 

.0000 

.0000 

.0000 

.0000 

.18)00 

.0000 

.188)0 

.0000 

1.0000 

.0000 

.0000 

-  .0001 

.01108 

.0081 

.1)118 

-  .0699 

—  .0681 

.0000 

.0000 

.0000 

.0000 

.18881 

.0188) 

.0000 

.18881 

.18881 

.0000 

.0000 

.18100 

.1888) 

.188)0 

.1327 

-  .18118 

.0008 

0.0000 

•fl 

r‘ 

sC 

II 

-  .0785 

-  .0883 

-  .0044 

-.1258 

1.1041) 

-.2381 

- .0255 

-  .11847 

.0532 

.2533 

.1003 

.7622 

.3554 

.8386 

-.8169 

1.1 42 II) +  1)5 

1.06901) +  05 

-  1.93420  +  04 

6.04630  +  04 

1.116571) +  05 

1.00881) +  05 

-  1.83880  +  04 

5.82620  +  04 

-  9.74I9D  +  05 

-4.44471) +  05 

1.65480  +  05 

-5.14870  +  05 

3.538  II) +  03 

2.70440  +  03 

-  3.32720  +  02 

1.47360  +  03 

-  1 .61011)  +  03 

-6.45451) +  02 

2.88230  +  02 

-  1.11250  +  03 

-  1.67271) +  04 

-  1.65310  +  04 

3.06050  +  03 

-9.518)80  +  03 

3.96521) +  01 

3.65220  +  01 

-7.38530  +  00 

2.32230  +  01 

4.8861 1) +  00 

6.13940  +  00 

-8.36400  -  01 

2.82940  +  00 

K  = 


-  .1)541 

- .0145 

-.1802 

-.1214 

.9174 

-  .5795 

-  1  9272 

4.18)79 

.5974 

-.1784 

.0396 

-  .4143 

-32.4650 

-  3.7504 

-  6.2656 

-  54.02M0 

.0628 

-.18)46 

.0335 

.3469 

-.6961 

4.1961 

6.2384 

-2.9799 

.2.346 

.0013 

-.1782 

24.1517 

-  18.3387 

-40.6498 

21.6354 

-  20.3557 

These  arc  as  follows. 

1  Rate  Command  (RC) 

2  Attitude  Command  w  ith  Attitude  Hold  (ACAH) 

3  Translation  Rate  Command  with  Position  Hold  (TRCPH) 

In  RC  systems,  attitude  must  diverge  from  trim  for  at  least 
4  seconds  following  a  step  input  command  In  ACAH.  a  con¬ 
stant  control  input  must  produce  a  proportional  angular  dis¬ 
placement  and  must  maintain  this  attitude  in  the  presence  of 
external  disturbances  In  TRCPH.  constant  control  input  must 
result  in  constant  translational  rate  and  the  rotorcraft  must  hold 
position  if  the  force  on  the  cockpit  controller  is  zero.  TRCPH 
systems  arc  preferred  in  nap  of  the  earth  maneuvers  in  fair  to 
poor  usable  cue  environments.  In  lair  usable  cue  environments 


considerable  concentration  is  required  for  the  pilot  to  perceive 
pitch  or  roll  attitude  and  lateral,  longitudinal,  or  vertical  trans¬ 
lation  rates  (Ret.  7).  The  use  of  cigcnstructurc  assignment  for 
the  design  of  RC  and  RCAH  has  been  discussed  in  a  previous 
paper  (Ref.  10);  therefore,  this  paper  will  concentrate  on  syn¬ 
thesis  of  a  TRCPH  using  eigenstructurc  assignment 

Both  clasMcal  (Ref.  1 1 )  and  modem  handling  qualities  lit¬ 
erature  (Refs.  7-9)  indicate  minimum  bandwidths  of  2  rad'see 
in  pitch,  roll,  and  yaw  rate  and/or  attitude  and  0.25  to  0.75 
rad/scc  in  vertical  velocity  arc  required  for  level  I  handling 
qualities.  In  addition,  coupling  between  lateral,  longitudinal, 
vertical  velocity,  and  vaw  modes  should  be  minimized  Band¬ 
width  requirements  of  0.2  to  I  rad  see  have  been  postulated 
for  TRCPH  (Ref.  12). 
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Fig.  I  (lain,  pole,  and  zero  locations  for  rotor  actuators  and  flexible 
modes:  l\i  for  real  eigenvalues.  uo»;  £i  for  complex  eigenvalues. 


Feedback  control  systems  must  maintain  stability  in  the  pres¬ 
ence  of  both  uncertainties  and  errors  in  the  mathematical  models 
used  lor  design  and  in  variations  in  system  parameters  durum 
actual  operations.  This  is  termed  stability  robustness.  Stability 
robustness  is  usually  specified  indirectly  in  terms  of  gam  and 
phase  margin.  Minimum  gain  margins  ol  6  db  and  minimum 
phase  margins  of  4?  are  ty  pical  stability  criteria  used  for  con¬ 
trol  lass  design. 


Control  Law  Synthesis 

The  control  lass  design  process  is  performed  in  tsso  stages 
tor  this  study  First,  a  full  state  regulator  is  developed  using 
eigenstrueture  assignment.  Even  though  good  decoupling  of 
the  closed  loop  eigenvectors  is  achieved,  control  coupling  ol 
the  helicopter  is  so  great  that  control  command  mixing  is  re¬ 
quired  Both  the  feedback  control  lass  and  the  mixer  are  de¬ 
signed  assuming  lull  state  feedback  The  feedback  control  lass 
cannot  be  implemented  directly  since  all  system  states  cannot 
be  measured  Thus  it  is  necessary  to  realize  the  control  lass  by 
means  ol  a  state  estimator  in  the  feedback  loop.  The  state 
estimator  is  synthesized  using  an  eigenstrueture  assignment 
technique  which  results  in  recovery  of  the  loop  transfer  prop¬ 
erties  of  the  lull  state  regulator. 
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Eig.  2  Open  limp  and  full  stale  class'd  limp  responses  to  a  side  velocity 
command: - open  loop. - full  stale  feedback. 


Full  Stale  KeRulalor 

The  leedbaek  control  law  is  a  linear  lunclion  of  the  state 

u  -  -  A'.i  i  3  i 

The  feedback  gain  matrix  K  is  selected  to  give  a  desired 
closed  loop  cigenvaluc'cigcnvcctor  configuration.  The  theory 
for  eigenstrueture  assignment  by  feedback  control  is  given  in 


Table  2  Open-Loop  Eigenvalues.  Eighth-Order  Model 


Open-Loop  Eigenvalues 

1.  -3.2610  +  .OOOOi  roll 

2.  -.9760  +  .OOOOi  pitch 

3.  .0820  ±  ,6296i  side  velocity 

4.  .1100  ±  ,SI47i  forward  velocity 

5.  -.2S88  ±  .0428i  yaw  vertical 


Open-lamp  Eigenvectors 


i. 

2. 

3. 

4. 

5. 

U 

-.0024 

.3605 

-  0710  *  .31031 

-.0462  +  ,4723i 

-  0244  -  .198.0 

\ 

-.0784 

.0102 

-.2171  +  .  12691 

-.0979  +  .  I297i 

.0616  +  ,0923i 

w 

-.0012 

.0525 

.0040  +  .024  li 

.0082  +  .034  li 

-.0949  +  .5468 

p 

- .8923 

.0373 

.1334  -  ,2352i 

- .0066  -  .I424i 

.0850  +  ,0084i 

q 

.0268 

-.6131 

-.0234  +  .28721 

.0666  +  ,2828i 

.0488  +  ,0270i 

r 

-.3407 

.3338 

.3382  +  ,4070i 

.4026  -f  ,2892i 

-  .7740  -  .  I489i 

.2844 

- .0767 

-.2699  -  .30261 

-.1982  -  . I093i 

-.0231  +  ,02.34i 

e 

- .0030 

.6104 

.4785  +  .07271 

.5862  -  .043 li 

-.0257  -  ,0?99i 
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Refs.  13-17.  The  details  of  the  application  of  the  theory  to  the 
problem  described  in  this  paper  is  given  in  Ref.  17:  therefore, 
this  section  will  be  limited  to  a  discussion  of  the  philosophy 
of  design  and  a  presentaton  of  the  results.  The  desired  eigen¬ 
values  and  eigenvectors  for  the  final  design  are  shown  in  Table  3. 
The  desired  eigenvalues  were  selected  to  satisfy  the  handling 
qualities  specifications  described  above.  The  roll  eigenvalue 
was  selected  to  give  a  roll  bandwidth  of  3.5  rad  sec.  a  pitch 
bandw  idth  of  2.9  rad  sec.  a  yaw  bandw  idth  of  3  rad  sec.  and 
a  vertical  velocity  bandw  idth  of  I  rad  sec.  The  eigenvalues 
associated  with  the  forward  and  side  velocities  were  selected 
to  be  complex  with  a  natural  frequency  of  0.9  rad  sec  and  a 
damping  factor  of  0.9.  This  resulted  in  bandw  idthx  of  about  I 
rad  sec  for  the  forward  and  side  velocities  as  well  as  the  vertical 
velocity. 

The  desired  eigenvectors  are  shown  in  Table  3  The  first 
desired  eigenvector  associated  with  roll  rate.  p.  is  made  equal 
to  unity  Since  the  roll  angle,  if.  is  the  integral  of  the  roll  rate, 
the  element  associated  with  roll  angle  is  the  inverse  of  the  roll 
eigenvalue  Also,  since  some  side  slip  is  inevitable  when  the 
helicopter  is  rolled,  an  arbitrary  non-zero  element  is  associated 
with  side  velocity  .  All  other  elements  of  this  eigenvector  arc 
zero  The  two  complex  eigenvectors  associated  with  the  side 
velocity  complex  roots.  -  0.802  -  0.388/.  have  side  velocity 
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elements  equal  to  one  and  arbitrary  non-zero  values  associated 
with  roll  rate  and  roll  angle.  All  other  elements  are  zero.  Thus 
the  lateral  modes,  roll  and  side  slip,  arc  decoupled  from  the 
longitudinal,  vertical  velocity  and  vaw  rate  modes  The  desired 
eigenvectors  associated  with  the  pitch  eigenvalue.  -2.9.  and 
the  forward  velocity  eigenvectors  associated  with  —0.801  — 
0.387/.  are  selected  in  a  similar  manner  to  decouple  the  lon¬ 
gitudinal  modes  from  the  lateral,  vertical  velocity  and  yaw  rate 
modes.  The  eigenvector  associated  with  the  vertical  velocity 
eigenvalue.  -  1 .0.  is  selected  so  that  all  components  are  zero 
except  the  vertical  velocity  ,  which  is  unity.  Similarly,  all  ele¬ 
ments  of  the  eigenvector  associated  with  the  vaw  rate  mode. 
-3.0.  arc  selected  as  zeros,  except  for  the  element  associated 
with  the  yaw  rate  which  is  chosen  to  be  unity. 

The  attainable  eigenvectors  for  a  unity  weighting  on  the  squared 
error  between  all  elements  of  the  desired  and  attainable  eigen¬ 
vectors  are  shown  in  Table  3.  Examination  ol  the  eigenvectors 
associated  with  roll  indicates  excellent  decoupling  between  this 
and  all  modes,  except  the  side  velocity.  The  pitch  mode  is  also 
decoupled  from  all  modes  except  the  forward  velocity  Yaw 
and  vertical  velocity  modes  are  also  decoupled  front  the  other 
modes  There  is  mild  yaw  coupling  in  the  side  velocity  eigen¬ 
vector  and  a  more  severe  vertical  velocity  coupling  in  the  tor- 
ward  velocity  eigenvector.  Responses  were  significantly 


Table  3  Closed  Loop  Desired  and  Attainable  Eigenvalues  and  Eigenvectors 


Closed-Loop  Desired  Eigentalues 

1.  -3.5000  +  .OOOOi  roll 

2.  -  2.9(881  +  .OOOOi  pitch 

3.  -  .8020  *  ,3880i  side  velocity 

4.  -  .8010  s  .3870i  forward  velocity 

5.  -1.18881  +  .OOOOi  heave 

6.  -3.1888)  +  .OOOOi  yaw 

Closed-Loop  Desired  Eigenvectors 

1.  2. 

3. 

4. 

5. 

6. 

u  .00(8) 

X 

.0000  + 

.OOOOi 

1.0000  + 

.OOOOi 

.1888) 

.0(881 

\  X 

.(MUM) 

1.00(H)  + 

.0000i 

.0000  + 

.OOOOi 

.1888) 

.(888) 

w  .0000 

.(HUM) 

.(HHM)  + 

.(MMK)i 

.0000  + 

.0000i 

1.0000 

.0000 

p  1 .0000 

.(MUM) 

X 

.0000  + 

.OOOOi 

•(H)(8) 

.0(M)() 

q  .(888) 

1.0000 

.(MMM)  + 

•OOOOi 

X 

.00(8) 

.00(8) 

r  .0(881 

.(MUM) 

.0000  + 

•OOOOi 

.0000  + 

.OOOOi 

.0(88) 

1.0(88) 

*  - .2857 

.0000 

X 

.0000  + 

.OOOOi 

.00(8) 

.(8881 

0  .0000 

-.3448 

.0000  + 

.0000i 

X 

.0(881 

.(MMM) 

****** 

UNITY  WEIGHTING 

*  *  *  * 

*  *  *  * 

*  *  *  * 

*  *  *  * 

********** 

********** 

***** 

***** 

Attainable  Eigenvectors 

u  .0024 

-.17  ' 

.0118  - 

,0307i 

.4254  - 

.OOOOi 

-.1003 

-.0112 

v  .0890 

-  .0235 

.4745  - 

.OOOOi 

.0131  + 

■0343i 

-.0181 

-  .0366 

w  0.0000 

- .0293 

-.0375  + 

•OI92i 

-  .2279  + 

•0045i 

.99.39 

-  .(8)28 

p  .9577 

.0041 

.3724  + 

.443  li 

.0389  + 

,0346i 

.0089 

-  .(8(65 

q  .0001 

.9357 

.0092  + 

.0564i 

-.4185  - 

,3999i 

-.0290 

-  .007.3 

r  -.0002 

-.00.36 

-.0282  + 

,0555i 

-  .018.1  - 

•0082i 

-.0014 

.9686 

d>  -  .2736 

.000.3 

-.5925  - 

,2726i 

-.0571  - 

.0171  i 

- .0089 

-.0.321 

0  0.0000 

- .3223 

-.0368  - 

.0559i 

.6197  + 

.  I999i 

.0290 

********** 

-.0142 

FINAL  WEIGHTING 

'  *  •  *  * 

*  •  *  * 

***** 

Attainable  Eigenvectors 


u 

.0024 

-.1381 

V 

.0890 

- .0235 

V* 

0.0000 

- .0293 

p 

.9577 

.0041 

q 

.0001 

.9357 

r 

-.0002 

-.0036 

6 

-.2736 

.0003 

e 

0.0000 

- .3223 

.0118  -  .030 Ji  .4143 

.4737  +  .OOOfli  .0093 

-  .0376  +  .0193i  -  .0024 

.3720  +  ,4473i  .0408 

.0084  +  .058 1  i  -.4350 

-.0028  +  ,0056i  -.0192 

-.5943  -  .27051  -.0606 

-.0369  -  ,0549i  .6440 


.OOOOi  -.1003  -.0112 

.03321  -.0181  -.0.366 

.0000i  9939  -  .0028 

.03771  .0089  -  .18165 

,4l54i  -  0290  -  .0073 

.00851  -.0014  .9986 

,OI94i  -.0089  -.0321 

,2075i  .0290  -  .0142 
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decoupled,  except  in  the  two  cases  mentioned  above.  This  was 
corrected  by  weighting  the  error  between  the  desired  and  at¬ 
tainable  vertical  velocity  element  in  the  forward  velocity  ei¬ 
genvector  by  a  factor  of  100.  The  yaw  coupling,  while  not 
serious,  was  also  reduced  bv  increasing  the  weighting  on  the 
error  in  the  yaw  direction  between  the  desired  and  attainable 
side  velocity  eigenvectors  by  a  factor  of  10.  The  attainable 
eigenvectors  for  this  set  of  weighting  terms  are  shown  in  Table 
3.  The  resulting  time  histories  for  lateral  step  input  for  this 
design  are  shown  in  Fig.  2  by  dashed  lines.  It  can  be  seen  that 
excellent  decoupling  has  been  attained  The  forward  velocity 
is  nearly  zero,  as  is  the  vertical  velocity.  Yaw  and  pitch  angular 
velocities  and  pitch  angle  are  also  nearly  zero.  The  roll  angle 
overshoots  its  steady  state  value  enough  to  produce  a  nearly 
constant  side  acceleration,  which  exists  until  the  commanded 
side  velocity  is  achieved.  The  roll  angle  then  decreases  until 
the  side  force  is  sufficient  to  counteract  the  steady  state  drag 
in  the  lateral  direction,  thereby  maintaining  constant  side  vc- 
locitx  Similar  results  were  obtained  for  longitudinal,  collec¬ 
tive.  and  directional  commands.  The  gain  matrix.  A.  for  this 
design  is  given  in  Table  I. 

In  smgle-inpul-single-output  (SISOl  systems,  stability  ro¬ 
bustness  is  measured  by  gain  and  phase  margins  obtained  from 
Bode  or  Nyquist  diagrams.  For  variations  in  the  input  gain  and 
phase,  the  MEMO  equivalent  ot  these  classical  stability  margins 
is  the  minimum  singular  value  (MSVi  of  the  return  difference 
matrix  |  /  *  A(.vlGl.vl)  for  s  =  tui  <  Rets  IK-21 1  The  Ml  MO 
gain  and  phase  margins  can  be  expressed  as 

Gain  Margin  (GMl  -  2(1  logi  I  tl  ~  MS V n 

Phase  Margin  I  PM  i  =  =  cos  ol  -  (MS\  F  2)  (4i 

For  a  minimum  singular  value  ol  unity,  the  gain  margins  arc 
infinity  and  6  db.  with  no  change  in  phase  at  the  input,  and 
the  phase  margins  are  r  60  deg  w  ith  no  change  in  gain  at  the 
input. 

The  minimum  singular  value  of  the  return  diftcrence  matrix 
tor  the  tull  state  design  is  plotted  (solid  linci  versus  frqueney 
in  Fig  4  The  minimum  singular  value  of  0.67  occurs  at  2  5 
rad  sec.  From  Eq  4  this  corresponds  to  gain  margins  of  4.63  db 
and  4  4?  db  and  phase  margins  of  =  34  deg.  Since  the  math 

model  of  the  helicopter  is  felt  to  be  reasonably  accurate  in  this 

frequency  range,  it  was  felt  that  these  margins  were  sufficiently 
large  to  guarantee  stability.  This  proved  not  to  be  the  case, 
however 

C  ontrol  Mixer  Design 

Even  with  the  modal  decoupling  in  the  closed  loop  system, 
the  control  coupling  resulted  in  excessively  coupled  responses 


Frequency  (rad/s) 

Fig.  3  Minimum  singular  valuts  vs.  frequency  for  full  state  feedback 
and  various  estimator  pole  locations. 


Frequency  (rad/s) 

Fig.  4  Minimum  singular  values  of  closed  loop  transfer  function  and 
maximum  singular  values  of  multiplicative  error  matrix  vs.  frequency . 


and  a  fixed  gain  control  mixer  was  required  to  convert  pilot 
inputs  to  control  inputs  ( note  the  control  mixer  described  below 
was  used  in  obtaining  the  results  in  Fig.  2l.  The  control  mixer 
was  designed  by  calculating  the  control  ti  required  to  produce 
a  given  steads  state  response  The  details  ol  this  design  arc 
given  in  Ref.  17.  The  resulting  control  mixer  gain  matrix.  //. 
is  given  in  Table  I .  The  closed  loop  system  equations  arc  now 

i  =  | A  -  «A|  i  *  BHb,  (5 1 

where  b,.  represents  pilot  command 

Compensator  Design 

The  control  law  developed  above  requires  knowledge  ol  the 
complete  state  vector  Since  only  measurements  of  roll,  pitch, 
and  yaw  rates  and  normal  acceleration  are  available,  j  stale 
estimator  is  required  in  the  feedback  loop.  The  control  is 

u  =  -  A  a  i  o  i 

where  X  is  the  estimate  of  the  state  given  by  the  state  estimator 

i  =  A. I  *  Bit  -*  /.( v  -  Ci  -  Dm  1 7 1 

It  is  well  known  that  an  estimator  such  as  given  by  Eq  7  may 
not  recover  the  stability  margins  of  the  lull  state  controller 
(Refs.  IK-21 ).  In  fact,  the  stability  margins  may  bo  very  poor, 
even  though  those  of  the  full  state  controller  may  be  excellent 
As  described  in  Ref  22.  an  estimator  which  recovers  the  lull 
state  loop  transfer  stability  properties  can  be  designed  by  ap¬ 
propriate  selection  ot  the  estimator  gain  matrix  L  such  that 
( I )  j  of  the  closed  loop  eigenvalues  ol  A  -  LC  approach  the 
finite  transmission  zeros  of  the  plant.  G't.v);  (2i  the  remaining 
n  -  j  closed  loop  eigenvalues  approach  infinity  ,  and  ( 3 1  the 
left  closed  loop  eigenvectors  of  A  -  LC  associated  with  the 
finite  eigenvalues  approaches  the  left  zero  direction  ol  the  finite 
transmission  zeros.  This  results  in  a  controller  in  which  the  lull 
slate  loop  transfer  matrix  Ft.vl  is  approximated  by  the  com¬ 
pensator  loop  transfer  matrix  Al.vlG't.vl  up  to  the  frequency  of 
the  infinity  poles. 

The  finite  transmission  zeros  of  G(.\)  were  calculated  av 
=  0  0 

e,  =  0  0 

=  -0.02147 
:4  =  -0  11224 
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In  designing  the  estimator,  the  gain  matrix  L  was  selected 
to  place  the  two  of  the  estimator  poles  at  :<  and  :4  and  two  of 
the  estimator  poles  at  —0.01  and  -0.012  (these  approximate 
the  transmission  zeros  at  zero).  Initially  the  remaining  estimator 
poles  were  placed  at  -  10.  -  12.  and  -5  ±  8.666:  (these 
approximate  the  transmission  zeros  at  infinity).  The  resulting 
singular  value  plot  is  shown  in  Fig.  4.  Stability  margins  arc 
almost  nonexistent.  The  'infinity”  pole  locations  were  in¬ 
creased  by  a  factor  of  ten  to  -  100.  -  120.  -50  ±  86.666/ 
The  singular  value  improves  substantially  but  is  still  unsatis¬ 
factory  at  high  frequencies.  Finally,  the  "infinity”  pole  lo¬ 
cations  were  increased  by  another  factor  of  eight  to  -  800. 
-960.  -400  ±  693.3/.  This  pole  configuration  resulted  in 
essentially  full  recovery  of  the  full  state  singular  values  over 
the  entire  frequency  range.  The  resulting  estimator  gain  matrix. 
L.  is  given  in  Table  3. 


Evaluation 

The  design  developed  above  using  the  eighth-order  model 
was  evaluated  using  the  37th  order  model  described  in  Fig.  I . 
The  MIMO  generalized  gain  and  phase  margins  of  9.63  and 
-  4  45  db  and  *  39  degrees  were  felt  to  be  adequate  to  provide 
stability  in  the  presence  of  unmodeled  dynamics  However,  the 
closed  loop  37th  order  system  exhibited  an  instability  near  30 
rad  s  due  to  coupling  between  the  main  rotor  collective: actuator 
dynamics  and  the  sensor  flexural  dynamics  as  measured  by  the 
normal  accelerometer  Instabilities  resulting  from  rotor  and  sen¬ 
sor  dynamics  have  been  noted  by  Chen  and  Hindson  (Ref.  23) 
and  Hall  and  Bryson  (Ref.  I ).  To  better  understand  the  source 
of  the  instability  the  multiplicative  error  matrix.  £(.v).  between 
the  8th  order  model  and  the  37th  order  model  was  computed 
at  various  frequencies.  This  matrix  is  defined  (Ref.  18)  as 

Gm(s)  =  (7*(.r)(/  4  £(.s)l  (8) 

To  be  assured  stability  in  the  face  of  modeling  errors  £(.v).  it 
is  known  (Ref  18)  that  at  all  frequencies 

<r(/  4  (A'*G«)'  'it?  5(£>  (9) 

A  plot  of  these  two  functions  is  shown  in  Fig.  5.  The  large 
error  peak  near  30  rad  s  crosses  the  minimum  singular  value 
curve  verifying  the  source  of  the  instability.  To  alleviate  this 
problem  one  might  increase  the  order  of  the  design  model  to 
include  dynamics  in  the  30  rad/s  range  and  completely  redesign 
the  compensator.  However,  a  filter  on  the  accelerometer  signal 
which  notched  out  the  set  of  complex  poles  at  30  rad's  and 
rolled  off  at  200  rad/s  was  implemented  and  was  found  to 
eliminate  the  instability. 

Transient  responses  to  pilot  commands  for  the  37th  order 
model  with  the  notch  filter  are  given  in  Fig.  5.  Figure  5  illus¬ 
trates  a  pure  side  velocity  maneuver  resulting  from  step  pilot 
lateral  input,  of  1  ft/sec.  In  this  case  longitudinal,  vertical 
velocity,  and  yaw  responses  are  minimal.  In  general,  the  re¬ 
sponses  in  Fig  6  arc  very  close  to  the  full  state  responses  of 
Fig  2  with  the  exception  of  small  lags  due  to  actuator  dynam¬ 
ics.  The  control  deflections  required  to  achieve  these  responses 
were  not  large  (Ref.  17). 

Closed  loop  transfer  functions  between  pilot  commands  and 
system  outputs  are  shown  in  Fig.  6.  These  figures  demonstrate 
(as  the  transient  responses  did)  that  the  closed  loop  system  is 
now  characterized  by  simple  decoupled  first  order  responses 
over  the  desired  bandwidth  range.  The  vertical  velocity,  for¬ 
ward.  and  sideslip  transfer  functions  arc  flat  until  about  0.8 
rad/sec  and  then  begin  to  roll  off.  The  yaw  transfer  function 
is  flat  until  about  I  rad/sec  then  rolls  off.  If  necessary  the 
bandwidth  on  these  transfer  functions  could  have  been  in¬ 
creased  by  increasing  the  natural  frequency  of  the  complex 
closed  loop  eigenvalues. 


•  il  ,5 

Time  (sec) 


*  c  ’  :  :■  j 

Time  (sec) 


Fig.  5  Closed  loop  lime  responses  to  a  side  velocity  command — 37th 
order  model. 


Conclusions 

The  eigenstructurc  design  techniques  described  in  this  paper 
provide  a  useful  method  for  the  design  of  control  laws  for 
helicopter  flight  control  systems  for  precision  hover  tasks.  Use 


Fig.  6  Bode  plots  of  closed  loop  transfer  functions  between  input 
commands  and  system  outputs. 
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of  MIMO  gain  and  phase  margins  based  on  minimum  singular 
values  of  the  loop  gain  did  not  predict  an  instability  due  to 
actuatonrotor  and  sensor  flexural  dynamics  in  the  collective 
control  loop.  This  instability  was  eliminated  by  filtering  the 
accelerometer  output  through  a  notch  filter.  The  use  of  an  error 
model  derived  from  an  estimate  of  the  dynamics  neglected  in 
the  design  model  did  predict  the  instability  and  it  is  recom¬ 
mended  that  this  approach  to  evaluating  stability  margins  be 
used  in  the  future.  The  state  estimator  in  the  feedback  loop 
required  high  gains.  This  might  cause  problems  in  actual  im¬ 
plementation  unless  the  sensor  outputs  are  filtered.  Also,  the 
effects  of  digital  implementation  of  the  control  law  s  might  cause 
difficulties  if  the  sampling  and  computational  rates  are  suffi¬ 
ciently  slow,  but  this  is  dependent  on  hardware  considerations 
and  is  not  addressed  in  this  paper 
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of  the  University  of  Minnesota. 
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